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Abstract 

We derive the semiclassical equations of motion of a transverse acoustical wave packet propa- 
gating in a phononic crystal subject to slowly varying perturbations. The formalism gives rise to 
Berry effect terms in the equations of motion, manifested as the Rytov polarization rotation law 
and the polarization-dependent Hall effect. We show that the formalism is also applicable to the 
case of non-periodic inhomogeneous media, yielding explicit expressions for the Berry effect terms. 
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I. INTRODUCTION 



The well-known analogy between the (linearized) equations of elasticity and Maxwell 
equations (see e.g. [lj) gives the advantage of predicting the acoustical analogue of many 
optical phenomena. In particular, polarization phenomena in optics find their counterparts 
in transverse acoustic waves. For example, the Rytov polarization rotation law [2] and the 
polarization dependent^ Hall effect in optics, which are manifestations of the Berry effect 



for photons [si, [4], [5], [f], Q, g[ ^, 10, [ll|, I12J, |l3], pertain to transverse acoustic waves as 
well ?], li, W^- The study of 'spin-orbit' interactions, that is the influence of intrinsic 
degrees of freedom (polarization/spin) on transport properties, has recently formed the 
subject of spin transport. Because of the dependence of transport on polarization/spin, 
particles in different polarization/spin states propagate along slightly different (semiclassical) 



trajectories. Exam 



effects in solids 



AH, 
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es of spin transport phenomena are the anomalous and the spin Hall 

inn n 

171 . Il8l . [19(, analogous effects for relativistic electrons 20], the optical 



13J and the acoustical [15[ polarization-dependent Hall effects. 
Recently, optical polarization transport in inhomogeneous media, in particular in_ pho- 
tonic crystals, has been studied via the semiclassical wave packet approximation |ll 
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211. 



The wave packet approximation had proved to be very successful in the study of electronic 



spin transport in solids, for which it was originally developed |16l. |22|. |23||. Here, via the same 
approach, we obtain analogous results for acoustical polarization transport in a phononic 
crystal (formed by the periodic variation of the density of the medium), subject to perturba- 
tions that vary slowly in space. Phononic crystals are novel materials that owe their name 
to the analogous effects of photonic crystals on light - they permit the passage of waves 
only at certain frequencies and can be, therefore, employed to control and confine sound 
|25| . l2q . l27l . |28| ] . A complication that arises in acoustical polarization transport in inho- 
mogeneous media is the coupling between the transverse and longitudinal waves (here caused 
by density variations), which means that polarization cannot be maintained. However, the 
coupling vanishes in the adiabatic limit of slow density variations, rendering polarization 
an adiabatic invariant. The formalism gives rise to a Berry phase term in the polarization 
evolution equation, describing the rotation of the polarization vector (the Rytov law), and 
a Berry curvature term in the equation of trajectory, deflecting the phonons depending on 
their polarization (the polarization dependent Hall effect). Such Berry effects have been pre- 
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viously established via a different semiclassical approach (the short-wave or post-geometric 
approximation) for transverse acoustical transport in non-periodic inhomogeneous media 



II. ACOUSTICAL BLOCH WAVES 

The dynamics of the displacement vector filed U(x, t) in an elastic medium is governed 
by (see e.g. [291 ) 

p%Ui = d 3 a l3 (1) 

(summation convention implied) where p(x) is the density of the medium and cr^- is the 
stress tensor. In an isotropic medium the latter takes the form 

Gij = X8ij d k U k + iiipiUj + djUi) 

A and n being the Lame coefficients (assumed constant throughout). Let us, thus, rewrite 
equation (CQ) as 

S£U = V(V ■ U) - ^ V x (V x U). (2) 

We decompose U into a solenoidal part, Uj_, and an irrotational part, Um, according to 
U = Uj_ + U|| (such a decomposition is, of course, always possible for an arbitrary vector). 
Since the volume change in a deformation is given by V • U, the solenoidal part of the 
displacement field does not involve any change in volume as it propagates in the body. It 
is thus termed the transverse displacement field, which has two independent components 
(polarization modes). The irrotational part, however, involves compressions and expansions 
in the body and is termed the longitudinal displacement field, having a single polarization 
mode. Thence, by taking the curl and the divergence of ([2]), we obtain 

V x (3 2 U ± - v\ V 2 U ± ) = -V x (<9 2 U,| - v\ V 2 U,|) = \Jv\ x V 2 Uj, 



V ■ (<9 2 U ± - v\ V 2 U ± ) = -V ■ (<9 2 U|| - v\ V 2 U,|) = -Vv\ ■ V 2 U ± (3) 



where 



p(x)' V P( x ) 

The curl and the divergence of the expressions in parentheses would vanish everywhere if 

the density were uniform, implying that the expressions must be zero identically, i.e., 

<9 2 U ± - v 2 ± V 2 V ± = 0, <9 2 U|| - w 2 V 2 U|| = (4) 



The transverse and longitudinal modes would, thus, propagate independently as waves with 
speeds v± and v», respectively. In presence of density variations, however, equations (j3J) 
are coupled. Nevertheless, introducing a slow (adiabatic) variation of the density, the right 
hand sides of may be neglected and the equations still decouple. This decoupling is 
vital in acoustical polarization transport, because otherwise, polarization cannot be main- 
tained (to be contrasted with optics). The adiabatic approximation eliminates the coupling 
between the transverse and longitudinal modes and renders polarization an adiabatic invari- 
ant. Therefore, in the adiabatic limit, (J4]) still applies where v± and vn are now local wave 
speeds, and we shall be interested only in the transverse wave equation, of course. 

The transverse index of refraction, n±, of an inhomogeneous medium is defined by v j_(x) = 
Cj_/nj_(x), where cj_ is the transverse wave speed in a reference homogeneous medium. For 
convenience, we shall take c± = 1. Introducing the monochromatic wave solution 



Ul(x,*)= e A C/ A (x,t) = 5> A « A (x)e-^ 

A=l,2 A 

e A being the (orthonormal) transverse polarization vectors, into the transverse wave equation 
yields the eigenvalue equation 
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n 2 , fx 



-V 2 u x = oo 2 u x (5) 



for polarized transverse waves. The differential operator in is not Hermitian. However, 
by introducing the Hermitian operator, 

<! = — 1 



71_l(x) n±(x) 

the equation can be written in the self-adjoint form 

QVa = (6) 

where 

V>a(x) = n_L(x)w A (x). 

Thus, with the weighting function nj_, eigenf unctions u x form a complete orthogonal set 
that is quasi (adiabatically) transverse. In a phononic crystal defined by imposing periodic 
condition on the refractive index through p(x), the eigenfunctions of (EI) are Bloch functions 
(Bloch's theorem), viz., 

V>An(x, q) = e iq - x w An (x, q) (7) 



where the corresponding eigenfrequencies, u; n (q), exhibit band structure as a result of the 
periodic condition. Here, q is the lattice momentum (restricted to the first Brillouin zone), 
n is the band index and u>a„(x, q) is periodic in x with the periodicity of the crystal lattice 
(viz. p). The acoustical Bloch waves (j7|) form a complete orthonormal set, 



and, thus, serve as a convenient basis for constructing phononic wave packets. In the follow- 
ing, we ignore the possibility of interband transitions and consider a single frequency band 
at a time, so that reference to the band index may be dropped. 

III. WAVE PACKET OF BLOCH PHONONS 

Let us proceed to the quantum level and introduce the ket notation, ip\(x, q) = (x\ip\(q)), 
w\(x, q) = (x.\w\(q)), etc. We consider a wave packet |\&) of Bloch phonons with a momen- 
tum k that is well defined on the scale of the Brillouin zone. Such a wave packet must be 
spread in real space over many primitive cells about a mean value r. In the semiclassical 
approximation, which describes the wave packet on a spatial scale much greater than its 
spread, the location of the packet may be specified by r. The semiclassical approach, thus, 
associates with each wave packet a position r and a momentum k (as well as a polarization 
state, see shortly), and predicts how these parameters evolve as the wave propagates in 
the medium. In the course of propagation, the wave packet may generally experience per- 
turbations in its neighborhood (such as deformation strain fields) that vary slowly within 
its width. Such slowly varying perturbations may be treated classically provided they are 
evaluated at the wave packet's mean position r. The semiclassical Hamiltonian H, which 
determines the dynamics of the wave packet via the time- dependent Schrodinger equation 
= H\ty) (h = 1), is then the sum of a periodic (quantum mechanical) Hamiltonian 
and the classically treated perturbations. Thus, H has the periodicity of the unperturbed 
Hamiltonian and involves r as a parameter. Whence r enters in the (Bloch) eigenstates and 
the corresponding eigenfrequencies parametrically. 

We, thus, construct the (normalized) wave packet according to 



dx V>L( X > q)^A'n'( X , q') = S X \'5 nn >5(ci ~ q') 




(8) 
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where g(q, t) and z\(t) are amplitudes with the normalizations 

/ dq\g(q,t)\ 2 = l 
Jbz 

£k(*)l 2 = i. 

A 

The momentum distribution \g\ 2 is assumed to be sharply peaked around its mean value k, 

k(t) = / dq \g\ 2 q 
Jbz 

which is regarded as the momentum of the wave packet. In the complex two-dimensional 
polarization space, the polarization state of the wave packet is represented by the parameters 
z\. The wave packet \^{t)) = \^(R(t))) is, thus, characterized by the parameters R = 
(R a ) = (r, k, z\). (For the sake of clarity, we reserve boldface letters for three-dimensional 
vectors, and middle alphabet indices for their components.) The location of the wave packet, 
which is specified by its mean position (\l/|x|\l/), is self-consistently given by 

r(t)= [ dqM 2 (V q S+(W^V q |W0), (9) 

where S(q, t) is the phase of the momentum amplitude, viz. g = \g\e~ lS , and |VT(r, q, z\)) = 
J2\ Z\\w\(r, q)) ((W^W 7 ) = 1)- In view of our assumption of a well localized wave packet in 
the momentum space, (jHJ) reduces to 

r(t) « V k S(k,t) + (W(R)\iV^\W(R)). (10) 

IV. DYNAMICS OF THE WAVE PACKET 

The dynamics of the wave packet is determined by the time-dependent Schrodinger equa- 
tion which is derivable from the Lagrangian 

£(|*),|*>,t)=i(*|*>-(*|^|*> 
via the Euler-Lagrange equation [30}. Thus, using (jSJ), 

C = f dq \g\ 2 (i(W\W) + d t S-u)^Yl z x {W{R)\id Zx \W{R)) + r • (W(R)\iV r \W(R)) 

JBZ A 

+d t S(k,t)-u(r,k) 

Since, 

d t S{k, t) = S(k, t)-k- V k 5(k, i)«5-k.(r- (W(R)\iV k \W(R))) 
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having used ([TO]) in the last expression, the Lagrangian becomes (up to an inconsequential 
total time derivative S) , 

£(R,R,t) « R-A(R) -k-r-u;(r,k), (11) 

where 

A(R) = {W(R)\iVz\W(R)) 
is the Berry connection (gauge potential) in the parameter space. It has components 

A n = (W(R)\id n \W(R)) = J2(MrMid n \w x ,(r,k)}z* x z x , 

AA/ 

A ki = (W{R)\id kt \W{R)) = J2(wx(r,k)\td ki \w Xf (r,k))zlz Xf 

AA/ 

A Zx = (W(R)\id Zx \W(R)) = iz$ 

The Berry curvature tensor (gauge field) associated with the Berry connection has compo- 
nents 

n ab {ii) = d Ra A h {R)-d Rb A a {R). 

Clearly, only A n and A ki give rise to non-vanishing components which are Vt nr . = (V r x 
A r ) fc = (^r)fc, tt k . k . = (V k x A k ) fc = (n k ) fc (i,j,k in cyclic permutations) and Q r . k .. By 
introducing matrix notation in the polarization space, these connections can be expressed 
more conveniently as 

A n = {z\A n \z), A ki = (z\A ki \z) 

where \z) is the column matrix with elements z x representing the polarization state ((z\z) = 
1), A r .(r, k) is the (Hermitian) matrix with elements (^(r, \t)\id ri \w x ,(r, k)) and likewise 
for A k .(r, k). The Lagrangian f[TT|) can be, thus, expressed as 

C ta i(z\z) + r • A r + k • A k - k- r - u. (12) 

Effective Lagrangians like f[T2T) . which determine the semiclassical dynamics of the wave 
packets, are central to the semiclassical wave packet approximation. The equations of motion 
that follow from (fl2|) are 

|i> =i(r- A r + k- A k )\z) (13) 
fi = (v g + k x n k )i - Q kir Tj (14) 
ki — (— V r u; + r x fl r )i + Q nkj kj (15) 
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where v s (r, k) = Vk^ is the group velocity of the wave packet. These are identical in form 
to the equations obtained in 21j for optical polarization transport in photonic crystals. 
Equation fll3p . which describes the time evolution of the polarization, has the solution 
\z(t)) = e^|z(0)), where 7 = J c (dr ■ A r + dk ■ A^) is the (Hermitian) Berry-phase matrix 
associated with the wave packet trajectory C in parameter space. The phase shift includes 
the contribution of the directional change of propagation (second term in the integrand) 
and thus generalizes the Rytov polarization rotation law. In the remaining two equations, 
the first terms on the right-hand side reproduce geometrical acoustics. The second term in 
( 1T4|) represents an anomalous velocity, which, being polarization dependent, gives rise to an 
acoustical polarization-dependent Hall effect, just as in the optical counterpart 12, 21]. 
The remaining terms owe to the existence of slowly varying perturbations. They are the 
second term in ffl5|) . which corresponds to a Lorentz force, and the third terms embodying 
the effect of mixed curvatures. 



V. NON-PERIODIC INHOMOGENEOUS MEDIA 

The above formalism is immediately applicable, just as in the optical counterpart, to the 

case of slightly inhomogeneous media with non-periodic refractive index, too. This, in the 

absence of perturbations, presents an example where the Berry connection/curvature can 

be calculated explicitly, reproducing the corresponding results of [ll|, Q for photons. In 

this case, a basis for the wave packet expansion is provided by |-0^(q)) = e* q ' x |e A (q)) so that 

\w\) — > |eA), eA being the polarization vectors. The latter can be written in the helicity 

basis as e A = -^(i e ± ii^), \ e (i^) being the zenithal (azimuthal) unit vector in the spherical 

coordinates of the momentum space. Thence, 

i /, \ 2 n cot^ 

A r (k) = 0, A k (k) = —jro-^ 



where (J3 is the Pauli matrix. This yields the Berry phase matrix 7 = cr 3 J c dk^^-, so 

cot t 



that right/left circularly polarized waves experience a geometric phase shift of ± Jcdk^ . . 
respectively. The non-vanishing curvature is given by 

Ok(k) = -^<*|*s|*) = -^(z(0)\a 3 \z(0)) 

where the final expression follows because 7 commutes with 03. Hence, Ok(k) = =Fp- for 
right/left circularly polarized waves respectively, which, according to fTbfl) . therefore propa- 



gate along different trajectories. These provide explicit examples of the Rytov law and the 
polarization-dependent Hall effect for transverse phonons. They coincide with the results of 
15( | obtained via the post-geometric approximation for non-periodic inhomogeneous media. 
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